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ABSTRACT 


Steele,  William  Gene,  M.S.A.A.E.,  Purdue  University,  December  1982. 
Aeroelastic  Stability  of  a  Forward  Swept  Wing  Aircraft  with  Wing  Tip  Stores. 
Major  Professor:  Terrence  A.  Weisshaar. 

When  studying  stability  of  an  aircraft,  a  quasi-steady  aerodynamic  model 
is  often  used.  In  this  study,  the  results  of  a  stability  analysis  of  an  unrestrained 
flexible  forward  swept  wing  aircraft  using  an  unsteady  aerodynamic 
formulation  will  be  compared  to  results  obtained  with  a  similar,  but  quasi¬ 
steady,  aerodynamic  formulation.  In  addition,  the  potential  effects  of  wing  tip 
stores  on  forward  swept  wing  aircraft  stability  will  be  explored. 
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ABSTRACT 


Steele,  William  Gene,  M.S.A.A.E.,  Purdue  University,  December  1982. 
Aeroelastic  Stability  of  a  Forward  Swept  Wing  Aircraft  with  Wing  Tip  Stores. 
Major  Professor:  Terrence  A.  Weisshaar. 

When  studying  stability  of  an  aircraft,  a  quasi-steady  aerodynamic  model 
is  often  used.  In  this  study,  the  results  of  a  stability  analysis  of  an  unrestrained 
flexible  forward  swept  wing  aircraft  using  an  unsteady  aerodynamic 
formulation  will  be  compared  to  results  obtained  with  a  similar,  but  quasi¬ 
steady,  aerodynamic  formulation.  In  addition,  the  potential  effects  of  wing  tip 
stores  on  forward  swept  wing  aircraft  stability  will  be  explored. 


INTRODUCTION 


Mission  requirements  for  modern  tactical  fighter  aircraft  dictate  the  need 
for  continued  research  in  advanced  fighter  technology.  Recent  studies  have 
shown  that  aeroelastic  divergence  of  swept  forward  wings  can  be  controlled  by 
the  application  of  aeroelastic  tailoring  of  composite  structures.  A  study  by 
N.J.  Krone  (1)  showed  that  tailored  composites  could  avoid  aeroelastic 
divergence  of  the  forward  swept  wing  with  little  or  no  weight  penalty.  The 
Defense  Advanced  Research  Projects  Agency  (DARPA)  funded  several  forward 
swept  wing  studies.  These  studies  identified  several  advantages  for  a  forward 
swept  wing  aircraft.  Krone  summarized  these  advantages  as  follows: 

1)  Configuration  flexibility 

2)  Significantly  higher  maneuver  -  L/D. 

3)  Lower  trim  drag-increased  supersonic  range  for  variable  sweep. 

4)  Lower  stall  speed. 

5)  Better  low  speed  handling. 

6)  Virtually  spin  proof. 

7)  Volume  benefits  -  lower  wave  drag. 

Also,  the  use  of  active  controls  for  flutter  suppression  and  enhanced  ride 
quality  has  become  more  feasible.  Thus  the  forward  swept  wing  (FSW)  fighter 
has  become  a  potentially  feasible  design  candidate  for  a  fighter  aircraft. 

This  study  has  two  principal  objectives.  First  it  will  examine  differences 
between  flutter  predictions  that  may  arise  from  the  use  of  quasi-steady  motion 
to  model  the  motion  dependent  aerodynamic  forces  on  the  flexible  FSW 
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aircraft  and  the  use  of  an  unsteady  representation  of  these  forces.  In  addition, 
the  study  will  identify  some  of  the  likely  effects  of  the  inclusion  of  wing  tip 
stores  on  the  stability  of  a  FSVV  aircraft. 

For  this  study  a  three  degree-of-freedom  model  is  used;  this  model  has 
body  freedoms  in  pitch  and  plunge  and  a  single  fundamental  bending  mode. 
Torsional  deformation  is  deleted  to  isolate  the  effect  of  bending  flexibility  on 
the  forward  swept  wing  stability  problem.  A  small  tip  missile  is  also  included 
in  the  analytical  model.  This  “store”  is  rigidly  attached  to  the  wing  tip. 
Important  parameters  to  be  varied  include:  wing  position  relative  to  the 
aircraft  c.g.;  and  tip  store  position. 
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BACKGROUND 


In  preliminary  stability  and  control  studies  the  analyst  focuses  on  the  rigid 
body  motions  of  the  aircraft.  In  contrast  to  this,  the  flutter  analyst  is 
concerned  primarily  with  the  elastic  vibratory  motions  of  the  aircraft  structure. 
For  a  rigorous  aeroelastic  study,  however,  the  rigid-body  and  the  elastic  modes 
should  be  considered. 

A  study  by  R.L.  Swaim  (2)  showed  that  there  is  a  great  deal  of 
interaction,  known  as  “Aerodynamic  Coupling”,  between  the  rigid-body  motion 
and  the  elastic  modes  of  a  flexible  vehicle.  This  interaction  can  drive  one  or 
more  of  the  modes  involved  unstable. 

In  a  study  by  Weisshaar  and  Zeiler  (3)  on  the  potential  vehicle  instability 
modes  of  a  flexible  FSW  aircraft,  quasi-steady  aerodynamics  were  used  to 
model  the  unsteady  aerodynamic  loads.  This  assumption  was  made  based  on 
the  fact  that  the  frequencies  of  vehicle  motion  for  this  type  of  analysis  were 
generally  quite  low. 

This  study  will  compare  the  aeroelastic  stability  results  of  a  FSW  aircraft 
using  both  unsteady  and  quasi-steady  aerodynamic  theories  to  model  the 
aerodynamic  loads.  Also,  it  will  investigate  the  effect  of  a  wing  tip  store,  and 
the  stores  chordwise  position  on  the  stability  of  the  FSW  aircraft. 

Future  fighter  aircraft  may  be  designed  with  forward  swept  wings, 
therefore  the  need  to  study  the  effects  of  wing  tip  stores  on  forward  swept  wing 
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METHODOLOGY 


Model  Parameters 

The  planform  geometry  of  the  aircraft  aeroelastic  stability  analysis  model 
is  shown  in  Figure  1.  Several  of  the  parameters  in  Figure  1  will  remain  fixed 
for  the  analyses  that  will  be  presented.  These  fixed  parameters  along  with  the 
mass  properties  of  the  aircraft,  are  given  in  Table  1. 

Table  1  Fixed  Aircraft  Parameters 


b  =  2.18  ft. 

m  —  1.656  slugs/ft.  (800  lb) 

/  =  15  ft 

Mt  =  496.89  slugs  (16,000  lb) 

d  =  4.5  ft 

n  -  0.11 

> 

II 

l 

CO 

o 

o 

Ia  =  138.67  lb-ft2/ft 

Figure  2  shows  a  typical  tip  missile  or  “store”  for  a  light-weight  fighter 
aircraft.  The  mass  properties  of  this  store  are  given  in  Table  2.  The  properties 
of  this  aircraft  model  were  considered  as  suitable  for  a  potential  light-weight 
fighter  aircraft.  Any  similarity  to  an  actual  aircraft  design  is  coincidental. 


Aeroelastic  Equations  of  Motion 

The  planform  geometry  of  the  aircraft  model  is  shown  in  Figure  1.  The 
flexible  wing  has  a  uniform  chord  and  sweep.  The  junction  between  the  wing 


Figure  2  Wing  Tip  Store 


Table  2  Store  Properties 

Ms  =  7.45  slugs  (240  lb) 

Is  =  170.13  lb-ft2 
0  =  0.30 

quarter  chord  line  and  the  side  of  the  fuselage  is  located  at  a  distance  x  aft  the 
aircraft  c.g.  The  swept  semi-span  of  the  wing  is  /.  The  model  has  three  degrees 
of  freedom.  These  are: 

1)  Upward  displacement  of  the  aircraft  c.g.,  denoted  as  w,  measured  positive 
upward.  This  is  the  plunge  degree  of  freedom. 
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2)  Pitching  of  the  fuselage  about  the  c.g.,  denoted  6,  measured  positive 
nose-up. 

3)  Bending  of  the  wings,  denoted  as  h,  measured  with  respect  to  the  wing 
root,  positive  upward. 

Aircraft  motion  is  restricted  to  be  symmetrical  with  respect  to  the  fuselage 
centerline. 

Three  equations  are  required  to  describe  the  symmetrical  motion  of  this 
aircraft.  One  equation  describes  the  vertical  motion  of  the  aircraft’s  center  of 
gravity.  A  second  equation  describes  the  motion  of  the  flexible  wing,  while  the 
third  equation  describes  the  pitch  motion  of  the  aircraft.  These  equations  were 
nondimensionalized  and  then  written  in  matrix  form.  The  form  of  these 
equations  is  as  follows: 

Qw 

Qn  (i) 

Qfl 

- 

The  elements  of  each  of  the  3x3  matricies  in  Eqn.  1  are  derived  in  Appendix  A. 
This  is  done  by  first  forming  expressions  for  the  aircraft’s  kinetic  energy,  strain 
energy  and  generalized  forces  due  to  the  aerodynamic  forces  acting  on  the 
airplane.  Lagrange’s  equations  are  then  used  to  generate  the  equations  of 
motion  of  the  airplane  when  it  is  perturbed  from  its  equilibrium  flight  path  by 
an  external  disturbance. 

The  elements  of  the  mass  matrix,  [M],  are  real  and  symmetrical.  Pitch 
and  plunge  are  not  coupled  dynamically  because  the  coordinates  chosen  to 
describe  aircraft  pitch/plunge  motion  are  located  at  the  aircraft  c.g.  However, 
bending  motion  is  inertially  coupled  together  with  both  pitch  and  plunge 


motion  in  the  mass  matrix.  Two  important  nondimensional  terms  appear  in 
the  mass  matrix  and  are  defined  as  follow's: 


The  elements  represent  the  wing  mass  to  fuselage  mass  ratio  and  store  mass  to 
wing  mass  ratios  respectively. 

The  elements  of  the  stiffness  matrix  [K],  are  derived  from  the  expression 
for  the  strain  energy  of  the  wing  in  terms  of  h,  the  tip  deflection,  and  an 
artificial  plunge  spring,  k. 

The  generalized  forces,  Qw,  Qh  and  Q?,  acting  on  the  aircraft  wings  and 
canard  arise  from  the  distributed  airloads  along  the  wings  and  a  single 
concentrated  force  on  the  canard.  These  loads  depend  on  wing  deformation 
and  aircraft  attitude  and  motion.  If  the  aircraft  is  given  an  arbitrary  virtual 
displacement,  consisting  of  <5w,  <5h,  60,  the  virtual  work  done  by  the  airloads  on 
the  wing  may  be  written  as  follows: 


6We  =  Qe60  +  Qh<5h  +  Qw<5w 


Using  Eqn.  4,  and  the  expressions  for  virtual  work  (see  Appendix  A)  the 
expressions  for  the  generalized  forces  in  terms  of  displacements,  w,  h  and  0  can 
be  formulated. 
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Methods  of  Solution 

The  solutions  of  the  time  dependent  equations  of  motion  are  of  the  form 
{£}est-  Making  this  substitution,  the  equations  of  motion  become: 

(s2(M)  +  [K]  +  q[A]]  {£}  =  {0}  (5) 

Eqn.  5  is  an  eigenvalue  problem  with  s  as  the  eigenvalue.  Matrix  [A]  is  the 
aerodynamic  influence  coefficient  (AIC)  matrix.  The  AIC’s  can  be  divided  into 
real  and  imaginary  parts.  They  are  functions  of  s  also. 

For  the  unsteady  aerodynamic  .presentation  the  AIC’s  are  functions  of 
reduced  frequency,  k.  However,  the  quasi-steady  representations  for  the  AIC’s 
are  independent  of  reduced  frequency. 

A  computer  program  was  developed  to  calculate  the  eigenvalues,  s,  for  the 
system  of  equations.  To  find  the  eigenvalues,  the  equations  of  motion  were  put 
into  the  following  form: 

!A]{X}  =  s{X}  (G) 

In  Eqn.  5,  s  is  the  eigenvalue  and  {X}  is  a  vector  of  displacements  and 
velocities  of  the  system.  To  convert  the  problem  from  that  given  in  Eqn.  5  to 
that  shown  in  Eqn.  6,  consider  the  following. 

The  equations  of  motion  for  this  system  can  be  written  as: 

|M|{q)  +  |RK<i>  +  |K|{q)  =  (0}  (7) 

A  new  variable,  is  now  introduced  to  reduce  this  second  order  system,  Eqn. 
7,  into  a  first  order  system. 


{q}  =  *{q} 
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Substitute  Eqn.  8  into  Eqn.  7.  The  resulting  equation  is: 

S[M]W  +  PITO  +  [K|{q>  =  {0}  (9) 

Rewriting  Eqn.  9  we  get: 

[Biro  -  [M]-‘  |K]{q>  =sro  (10) 


Eqn.  8  can  be  rewritten  as: 


[I]  ft}  =  sq 


(11) 


Combining  Eqns.  10  and  11,  the  equations  of  motion  can  be  put  into  the  form 
shown  in  Eqn.  6.  In  this  case: 


(A]  = 


(0]  [1] 

-[M)_,(K]  — [M]~ 1  [B] 


(12) 


and 


{X}  = 


q 

q 


(13) 


Matrix  [A]  is  a  6x6  matrix.  There  are  six  roots  to  this  equation.  They  may  be 
either  complex  conjugates  or  real.  The  complex  conjugate  eigenvalues  are  of 
the  form: 


s  —  (14) 

Where  a  represents  the  system  damping  and  u>  is  the  frequency  of  oscillation. 

These  roots  may  be  plotted  in  a  “root  locus”  format  using  <r  and  ju>  as 
coordinates.  The  roots  with  negative  frequencies  are  omitted  from  these  root 
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locus  plots  since  they  are  merely  reflections  of  the  other  complex  conjugate 
roots.  However  the  roots  with  a  zero  frequency  are  of  interest.  The  root  locus 
plot  is  a  valuable  tool  for  determining  system  stability.  Thus  stability  can  be 
studied  by  examining  the  behavior  of  roots  in  the  “s-plane”.  Instabilities  occur 
for  values  of  damping,  <r,  greater  than  or  equal  to  zero. 

Two  solution  techniques  were  used  to  solve  the  flutter  problem.  The  first, 
using  an  unsteady  representation  of  the  aerodynamic  forces,  is  called  the  P-K 
method.  The  other  method  is  a  variation  of  P-K  method  using  quasi-steady 
aerodynamics.  A  detailed  description  of  these  methods  is  given  in  Appendix  B. 

Several  aircraft  parameters  are  to  be  varied.  These  parameters  are  the 
wing  position,  5T,  wing  tip  store  position,  5TS,  and  the  store  to  wing  mass  ratio, 
/?.  Several  different  situations  were  considered  and  are  listed  in  Table  3.  To 
begin  with,  the  stability  of  the  aircraft  without  tip  stores  was  studied.  These 
so-called  “clean  wing”  cases  are  compared  to  the  cases  with  stores  to  determine 
the  effects  the  addition  of  the  wing  tip  stores  have  on  the  stability  of  the 
system.  Also,  calculations  were  made  using  quasi-steadv  and  unsteady 
aerodynamics  for  identical  geometrical  and  inertial  input  parameters.  A 
comparison  of  the  results  obtained  from  these  aerodynamic  formulations  will  be 
made  to  determine  differences  in  the  calculated  flutter  speed  of  the  aircraft. 


Table  3  Clean  Wing  and  Wing/Store  Cases 


Case 

No. 

P 

Store 

Loc. 

Quasi- 

Steady 

Unsteady 

1 

m 

X 

mi 

H 

m 

X 

B 

0.45 

X 

Clean 

5 

0.40 

X 

Wing 

6 

0.40 

X 

7 

0.35 

X 

8 

0.35 

X 

9 

0.30 

X 

10 

0.30 

X 

11-15 

0.50 

0.30 

0-c 

X 

Wing 

16-20 

0.50 

0.30 

0-c 

X 

with 

21-25 

0.40 

0.30 

0-c 

X 

Store 

26-30 

0.40 

0.30 

0-c 

X 

31-35 

0.30 

0.30 

0-c 

X 

36-40 

0.30 

0.30 

0-c 

X 

41 

0.30 

0.30 

0 

X 

42 

0.30 

0.60 

0 

X 

♦Stores  were  varied  from  leading  edge,  0,  to  trailing  edge,  c, 
at  1/4  c  increments. 


RESULTS  AND  DISCUSSION 


The  first  case  we  will  consider  is  an  aerodynamically  stable 
wing/body/canard  combination  without  tip  stores.  This  airplane  has  a  30° 
forward  swept  wing,  ft  =0.11,  and  a  cantilever  wing  flexural  frequency 
cj0  =  68rad/sec.  The  wing  root  fuselage  junction  is  aft  of  the  aircraft  c.g.  at 
the  position  X  =  0.35.  (see  Figure  1).  Note  that  the  barred  quantities,  such  as 
T,  refer  to  a  parameter  that  is  divided  by  the  wing  semi  span,  /. 

Figure  3  shows  a  root  locus  plot  of  the  eigenvalues  for  a  quasi-steady 
aerodynamic  analysis  of  this  configuration  as  airspeed,  V,  increases.  Each  point 
on  the  plot  corresponds  to  a  different  airspeed.  These  values  for  airspeed  are 
shown  nondimensionalized  with  respect  to  a  reference  speed,  VDC,  the  airspeed 
at  which  the  clamped  wing  encounters  aeroelastic  divergence. 

Two  curves  are  plotted  in  Figure  3.  The  first  represents  the  behavior  of 
the  pitch  mode  (also  called  the  short  period  mode )  as  V/VDC  increases.  The 
frequency  of  the  pitch  mode  first  increases  while  pitch  damping  also  increases. 
However,  near  V/VDC  =  0.90,  frequency  begins  to  decrease  until  the  root  locus 
plot  intersects  the  real  axis.  At  this  point,  the  root  divides  and  moves  along 
the  real  axis. 

The  coupled  bending  mode  decreases  in  frequency  and  increases  in 
damping,  <r,  until  a  velocity  near  V/VDC  =  0.80.  The  bending  root  then 
moves  toward  the  right  half  plane  as  airspeed  increases.  At  VF/VDC  =  0.949, 
the  bending  mode  becomes  unstable.  Aeroelastic  coupling,  involving  the  two 


Root 


(l  /sec) 


5cus  Plot  for  3T  —  0.35,  Quasi-Steady  Aerodynamic 


16 


- 


1 


t; 


rigid-body  freedoms  and  the  single  bending  mode,  leads  to  a  dynamic 
instability  of  the  airplane  referred  to  as  body-freedom  flutter  (3).  Although  the 
traditional  scheme  for  identifying  aeroelastic  modes  (bending,  pitch)  has  been 
used  here,  these  are  not  normal  modes.  In  fact  they  each  contain  large 
amounts  of  bending,  plunge  and  pitch. 

One  other  root  is  present  in  the  analysis,  but  is  not  shown.  This  root 
involves  the  plunge  mode  and  is  located  at  and  near  the  origin  along  the  er  axis. 
Since  this  root  is  of  no  significance  to  the  aeroelastic  stability  analysis,  it  is 
omitted  from  the  root  locus  plots  to  follow. 

Figure  4  is  a  root  locus  plot  for  the  same  aircraft  configuration.  However 
unsteady  aerodynamic  theory  is  used  to  model  the  aerodynamic  forces.  In  this 
figure  the  coupled  bending  mode  is  well  damped  at  all  airspeeds.  The  coupled 
pitch  mode  frequency,  w,  and  damping,  er,  increase  with  increasing  airspeed 
until  a  velocity  slightly  above  V/VDC  =  0.85.  The  pitch  root  then  moves 
toward  the  right  half  of  the  s-plane  as  airspeed  increases.  At  VF/VDC  =  0.96 
the  pitch  mode  becomes  unstable  at  a  frequency  of  18.21  rad/sec  and  a  reduced 
frequency  k  =  0.0203. 

The  most  noteworthy  difference  between  these  two  root  locus  plots, 
Figures  3  and  4,  is  the  apparent  change  in  the  origin  of  the  flutter  mode.  The 
the  flutter  speed,  as  calculated  from  the  quasi-steady  representation,  is  slightly 
lower  than  the  flutter  speed  calculated  with  the  unsteady  aerodynamic 
representation,  (see  Figure  4).  However,  although  the  ancestry  of  the  flutter 
mode  is  different,  the  coupled  instability  mode  is  nearly  the  same  in  both  cases. 

The  next  case  considered  is  an  unstable  wing/body  combination.  This 
airplane  is  identical  to  the  one  given  in  the  previous  discussion,  except  that  the 
wing  root  position  is  moved  closer  to  the  aircraft  c.g.  to  a  position  5T  =  0.30 
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As  before  two  root-locus  plots  are  of  interest. 

Figure  5  is  a  root  locus  plot  for  this  configuration  using  a  quasi-steady 
representation  for  the  aerodynamic  loads.  The  coupled  wing  bending  mode 
damping  increases  while  frequency  decreases  with  airspeed  until  about 
V/VDC  =  0.80.  The  bending  root  then  moves  toward  the  right  half  plane  as 
airspeed  increases.  At  Vp/VDC  =  1.003  the  wing  bending  mode  becomes 
unstable.  For  the  coupled  pitch  mode  the  roots  are  non-oseillatory.  Some  of 
these  roots  pass  into  the  right  half  plane.  This  indicates  a  “static”  or  aperiodic 
instability  in  the  pitch  mode.  To  further  illustrate  this  instability,  Figure  6 
shows  damping,  <r,  vs.  non-dimensional  velocity,  V/^DC  f°r  e^ch  mode.  This 
figure,  like  the  root-locus,  indicates  an  instability  when  values  for  a  are  greater 
than  zero.  As  seen  in  Figure  6,  there  are  two  coupled  pitch  modes,  one  stable, 
the  other  unstable.  This  figure  also  shows  the  bending  mode  root  crossing  the 
zero  damping  line,  indicating  the  onset  of  flutter,  and  the  unstable  pitch  root 
re-stabilizing  in  the  same  vicinity. 

Figure  7  is  a  root  locus  plot  for  this  same  aircraft  configuration  using 
classical  Theodorsen  unsteady  aerodynamic  theory  to  model  the  aerodynamic 
forces.  As  before,  the  coupled  wing  bending  root  damping  increases  and 
frequency  decreases  with  airspeed  until  near  V/VDC  =  0.90.  This  root  then 
moves  toward  the  right  half  plane  as  airspeed  increases.  At  Vp/VDC  =  1.033 
the  wing  bending  root  becomes  unstable.  This  flutter  speed  is  slightly  greater 
than  that  obtained  from  using  quasi-steady  aerodynamics.  The  most 
significant  difference  between  Figures  5  and  6  is  seen  in  the  pitch  mode.  In 
Figure  5  the  pitch  mode  was  non-oscillatory  and  statically  unstable  at  a  low 
airspeed.  However,  when  using  unsteady  aerodynamics,  the  pitch  mode 
becomes  oscillatory  and  the  static  instability  disappears. 
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Modal  Danping  vs.  Velocity  Ratio  when  3T  =  0.30,  Quasi- 
Steady  Aerodynamic  Theory 


Another  indication  of  differences  resulting  from  the  quasi-steady  and 
unsteady  aerodynamic  formulations  is  seen  by  examining  a  plot  of  the  aircraft 
flutter  speed  as  the  wing  position,  JT,  changes.  As  shown  in  Figure  8,  the  flutter 
speeds  as  predicted  by  quasi-steady  and  the  unsteady  aerodynamic  theory 
intersect  each  other  near  3T  =  0.37.  This  indicates  that  the  quasi-steady 
aerodynamic  formulation  does  not  always  give  a  larger  value  for  the  flutter 
speed. 

When  stores  are  added  to  the  wing  tips  of  the  aircraft,  very  little  change 
in  the  root  locus  plots  was  observed  as  the  store  c.g.  location,  5TS,  was  moved 
fore  and  aft  along  the  chord  of  the  wing.  Figure  9  is  a  root  locus  plot  for  a 
wing  position,  S',  of  0.40,  store  to  wing  mass  ratio,  /?,  of  0.30  and  the  store 
location  on  the  wing  tip,  3TS,  of  0.0,  the  mid-chord  position.  Quasi-steady 
aerodynamics  were  used  to  model  the  aerodynamic  forces.  Flutter  occurs  in 
the  coupled  wing  bending  mode. 

Figure  10  is  a  root  locus  plot  for  the  identical  aircraft  using  an  unsteady 
aerodynamic  representation  of  the  aerodynamic  forces.  In  this  figure,  flutter 
occurs  in  the  coupled  pitch  mode.  Figures  11,  12,  and  13  arc  plots  of  flutter 
speed  vs.  store  position.  In  each  of  these  figures  the  quasi-steady  flutter  speed 
decreases  as  the  store  is  moved  from  leading  edge  (L.E.)  to  trailing  edge 
(T.E.).  In  Figures  11  and  12  the  flutter  speed  found  with  the  unsteady 
aerodynamic  representation  increases  as  the  store  is  moved  from  L.E.  to  T.E. 
This  trend  is  opposite  to  that  seen  with  the  quasi-steady  aerodynamic- 
representation.  Also,  as  was  seen  for  the  clean  wing  cases,  quasi-steadv 
aerodynamics  docs  not  always  give  a  lesser  value  for  the  flutter  speed.  In 
Figure  11,  the  quasi-steady  flutter  speed  is  higher  than  the  unsteady  flutter 
speed,  but  in  Figures  12  and  13  the  opposite  is  true. 
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Figure  9  Root  Locus  Plot  for  5T  =  0.40,  Quasi-Steady  Aerodynamic 

Theory,  with  Store  c.g.  at  Wing  Midchord  (JTS  =  0) 
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Another  interesting  difference  between  the  results  obtained  with  the  two 
aerodynamic  formulations  is  found  from  a  comparison  of  the  clean  wing  flutter 
speeds  with  the  flutter  speeds  obtained  with  wing  tip  stores.  Figures  14  and  15 
are  plots  of  flutter  speed,  Vp/VDC  vs.  wing  position,  7,  for  the  clean  wing  cases 
and  for  the  cases  with  a  store  located  on  the  midchord.  In  these  figures,  it  is 
evident  that  the  flutter  speed  decreases  when  stores  are  added  to  a  clean  wing. 
Figure  14,  for  quasi-steady  aerodynamics,  shows  that,  as  the  wing  moves  aft  of 
the  aircraft  c.g.,  the  differences  between  the  flutter  speeds  of  the  clean  wing 
and  the  wing  with  a  tip  store  increase.  However,  with  unsteady  aerodynamics, 
Figure  15  shows  that  the  difference  between  the  flutter  speeds  remains 
relatively  constant. 

It  was  noted  earlier  in  this  study  that  the  use  of  different  aerodynamic 
theories  to  represent  the  aerodynamic  forces  of  the  aircraft  model  changed  the 
flutter  mode.  This  flutter  mode  change  also  occurs  when  wing  tip  stores  are 
added.  Figure  9  is  a  root  locus  plot  for  3T  =  0.40,  (3  —  0.30  and  3TS  =  0. 
Flutter,  for  this  aircraft  geometry,  occurred  in  the  coupled  wing  bending  mode. 
Figure  16  is  a  root-locus  plot  for  the  same  aircraft  configuration,  with  the 
store  removed.  Flutter  now  occurs  in  the  coupled  pitch  mode. 

For  a  forward  swept  wing,  a  situation  could  arise  such  that  the  wing 
center  of  pressure  lies  ahead  of  the  aircraft  c.g.  For  an  aircraft  with  a  canard, 
this  configuration  will  be  statically  unstable.  Figure  5  is  an  example  of  a.  root 
locus  plot  for  such  a  configuration.  Figure  17  is  a  root  locus  plot  for  the  same 
aircraft  configuration  but  with  the  addition  of  a  tip  store  on  the  leading  edge  of 
the  wing.  The  aircraft  parameters  are  /?  =  0.30,  and  x  =  0.30.  Quasi-steady 
aerodynamics  were  used  for  these  results.  Note  that  with  the  addition  of  the 
store,  this  configuration  is  no  longer  statically  unstable  because  the  aircraft  c.g. 
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Figure  14 


Flutter  Speed  Ratio  vs.  Wing  Position  with  and  without 
Stores  Placed  at  the  Midchord,  Quasi-Steady  Aerodynamic 


Figure  15 


Flutter  Speed  Ratio  vs.  Wing  Position  with  and  without 
Stores  at  Midchord,  Unsteady  Aerodynamic  Theory 


has  been  moved  forward  by  the  addition  of  the  store  on  the  wing  tip.  In  fact, 
the  pitch  mode  went  from  non-oscillatory  to  slightly  oscillatory.  Also  the 
flutter  speed  was  reduced  from  VF/VDC  =  1.003  without  a  store  to 
VF/VDc  =  0.98872  with  the  store.  To  further  illustrate  this  finding,  another 
root  locus  plot  was  generated  for  the  same  aircraft  parameters  as  the  previous 
case.  However,  the  store/wing  mass  ratio,  /?  has  been  doubled.  Figure  18  is  a 
root  locus  plot  for  x  =  0.30  and  /?  =  0.60.  Again  quasi-steady  aerodynamics 
were  used.  Note  the  increase  in  frequency  for  the  pitch  mode.  Again,  as 
before,  the  addition  of  the  wing  tip  store  has  further  reduced  the  flutter  speed 
from  Vp/VDC  =  0.989  to  Vp/VDc  =  0.978.  With  the  addition  of  the  store  on 
the  wing  tip,  the  aircraft  c.g.  moved  sufficiently  forward  ahead  of  the  wing 
center  of  pressure  to  eliminate  the  static  instability. 
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Figure  18 


Root  Locus  Plot  for  X  —  0.30,  /?  =  0.60,  Quasi-Steady 
Aerodynamic  Theory  with  Store  c.g.  at  Leading  Edge 


SUMMARY  AND  CONCLUSIONS 


One  of  the  purposes  of  this  study  was  to  identify  potentia1  differences  in 
the  stability  characteristics  of  an  idealized  forward  swept  wing  aircraft  model 
when  quasi-steady  and  unsteady  aerodynamic  formulations  are  used.  Also,  the 
effect  of  the  inclusion  of  wing  tip  stores  on  the  stability  of  the  forward  swept 
wing  aircraft  model  was  considered.  Because  of  the  simplicity  of  the  model  and 
the  limited  number  of  cases  analyzed,  general  conclusions  relevant  to  all 
forward  swept  wing  aircraft  cannot  be  made.  However,  from  these  preliminary 
studies,  some  conclusions  can  be  drawn. 

Several  differences  were  noted  in  the  stability  characteristics  of  the  aircraft 
model,  depending  upon  whether  or  not  quasi-steady  or  unsteady  aerodynamic 
theory  was  used.  Flutter  occurred  in  a  different  mode  depending  on  the 
aerodynamic  formulation  used.  A  statically  unstable  aircraft  configuration 
using  quasi-steady  aerodynamics  was  found  to  be  stable  when  unsteady 
aerodynamics  were  used.  It  was  also  shown  that  quasi-steady  aerodynamics 
does  not  always  give  lower  values  for  the  flutter  speed  than  does  the  unsteady 
aerodynamic  representation. 

When  an  aeroelastic  stability  analysis  is  performed  on  an  aircraft,  it  is 
necessary  to  consider  all  the  aircraft  modes  to  determine  preciseh  which  mode 
caused  the  instability.  Also,  because  of  the  differences  in  the  stability 
characteristics  encountered  when  using  the  quasi-steady  and  unsteady 
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formulations  for  the  aerodynamic  loads,  an  accurate  representation  of  these 
loads  is  essential. 

When  stores  were  added  to  the  wing  tip,  the  flutter  speed  generally 
decreased  for  quasi-steady  aerodynamics  as  the  store  was  moved  from  leading 
edge  to  trailing  edge.  However,  for  one  of  the  wing  positions  studied  this  trend 
was  reversed  when  unsteady  aerodynamics  were  used.  It  was  also  found  that, 
with  quasi-steady  aerodynamics,  the  differences  between  the  flutter  speeds  for 
the  clean  wing  and  the  wing  with  stores  increased  as  the  wing  was  moved  aft  of 
the  A/C  c.g.,  while,  with  unsteady  aerodynamics  the  differences  between  the 
flutter  speeds  remained  relatively  constant.  Also,  when  stores  are  added  to  the 
wing  tip  of  a  statically  unstable  FSW  aircraft  geometry,  it  was  found  that  this 
instability  disappeared.  Generally,  however,  it  was  found  that  the  flutter  speed 
was  lower  when  stores  were  added  to  the  wing  tips  regardless  of  the 
aerodynamic  representation  used  for  this  study. 
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APPENDIX  A 

Development  of  Perturbed  Equilibrium  Equations 
With  the  Addition  of  Wing  Tip  Stores 


This  appendix  presents  the  development  of  the  equations  of  motion  for  a 
three  degree-of-freedom  forward  swept  wing  model  with  wing  tip  stores.  A 
planform  view  of  this  model  is  shown  in  Figure  A-l.  The  degrees-of-freedom 
are:  w,  aircraft  plunge;  6,  aircraft  pitch  attitude;  and,  h,  wing  bending  or 
flexure. 


Figure  A-l  Planform  Geometry  and  Nomenclature 
The  equations  of  motion  are  developed  by  using  Lagrange’s  equations, 
together  with  an  assumed  deflection  shape  for  wing  bending  vibration. 
Expressions  for  the  kinetic  energy,  T,  the  strain  energy,  U,  and  expressions  for 


I 
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the  generalized  aerodynamic  forces  associated  with  each  degree  of  freedom,  will 
be  formulated. 

The  kinetic  energy  of  the  aircraft  and  stores  is  written  as: 


T  =  7  |Mf  (w  +  e0)2  +  IF  0s] 

mt 


(A-l) 


+  2(— )f  f  p(w  -  (z  +  y  sinA  +  £  cosA)0  +  h]2dfdy 


+  7  (2MS)  h,2  +  |  (2IS)  S2 


The  first  bracketed  term  Eqn.  A-l  represents  the  kinetic  energy  of  the  fuselage; 
the  integral  term  accounts  for  the  kinetic  energy  of  the  mass  distributed  over 
the  wing  surface  area  S.  The  last  two  terms  represent  the  kinetic  energy  of  the 
wing  tip  stores.  The  coordinates  (y,£)  define  the  location  of  a  mass  point  on 
the  wing  with  respect  to  the  wing  root  origin  located  a  distance  z,  aft  of  the 
aircraft  c.g.  The  vertical  velocity  of  the  store  c.g.,  hs,  is  given  by  the 


expression: 


hs  =  w  +  h  -  (z  +  /  sin  A  +  xs)0 


(A-2) 


This  analysis  will  be  restricted  to  a  uniform  planforrn  wing  with  constant 
mass  per  unit  area  distributed  uniformly  about  the  y-axis.  With  this 
restriction  Eqn.  A-l  becomes: 


T  =  w2  +  2  MFe  w 0  +  (MFe2  +  IF)^J  +  m £  [2w2  +  h2 


+  ^(z2  +  y2sin2A  +  r2  cos2A)]dy  +  m  /  [2  wh  -  2  wt?(z  +  y  sinA) 
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-  2h 9  (z  +  y  sinA)  +  (2zy  sinA)  0s]  dy  +  [w2  +  2hw 

+  h2  -  2(z  +  /  sinA  +  xg)  #w  -  2  (z  +  /  sinA  +  xg)h# 

+  (z  +  /  sinA  +  Xg)2#2  4-  Is  (P  (A-3) 

In  Eqn  A-3,  use  has  been  made  of  the  fact  that  J £d£  is  to  be  zero  for  a 
uniformly  distributed  mass.  In  addition,  the  moment  of  inertia  per  unit  length 
of  the  wing  sections,  about  the  y-axis,  is  IQ  =  mr2  and  the  moment  of  inertia 
of  the  store  is  Is  =  Ms  r82 

The  mass  matrix  for  the  model  is  developed  by  first  assuming  a  mode 
shape  wing  flexure,  h(y,  t).  This  shape  function,  f(y),  is  given  as: 

f(y)  =  J  [6(y /If  +  4(y /If  +  (y/o'l  (A-4) 

so  that 

h(y,t)  =  E(t)f(y)  (A- 5) 

the  kinetic  energy  expression  is  now  given  as: 

T  =  “{(MF  +  2  ml)w2)  +  ±{IF  +  MFe2} 

+  2  m/(z2  +  sin2A  +  r2  cos2A  +  zl  sinA)}#2 

o 

+  -^{2mf  f2dy}  E2  +  {2m/  fdy}wE  (A-6) 

i  / 

+  (M^e  -  2m/  (z  +  —  l  sinA))#w  -  {2m /  (fz  4-  fy  sinA)  dy)#E 
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+  (z  +  /  sinA  +  x,)2#2]  +  M^,2#2 


Since  the  airplane  motion  is  described  by  the  coordinates  w,  9,  and  E,  the 
mass  matrix  for  the  equations  of  motion  for  the  model  is  developed  by  using 
the  standard  definition  of  Mjj  in  terms  of  partial  derivatives  of  T. 


—  ~~~~  —  Mp  +  2m/  +  2  Mg  —  M-p  +  2Mg 

dv." 


d2T  J 

M12  =  JLrL^  =  2m/  fdy  +  2  Ms 

dwdE  0 


M13  —  -  d  — =  Mpe  -  2  m/(z  +  —  /  sinA) 
dwdO  2 


-  2MS  (z  +  /  sinA  +  xs) 


(A-7) 


(A-8) 


(A-9) 


Because  the  reference  system  is  centered  at  the  aircraft  center  of  mass,  Mj3  is 
equal  to  zero.  From  statics  the  following  equation  for  e  is  found: 


e  —  2  77-  (z  +  ~  /  sinA)  +  277—  (z  +  /  sinA  +  xs)  (A-10) 
Mp  2  Mp 


The  remaining  mass  matrix  terms  are: 


Moo  =  d-r-T  =  2  m  /  f2  dy  +  2  Ms 
dh~  0 


(A-ll) 


Mo3  =  ~  =  -  2m /  f(z  +  y  sinA)  dy 

dKd9 


(A-12) 


-  2  Ms  (z  +  /  sinA  +  xs) 

M33  =  -~o-  =  If  ■*"  Mpe2  +  2  ml  (z2  +  x/  sinA  +  —  /2  sin2A  +  r2  cos2  A) 

ar  3 
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+  2  Ms(z  +  /  sinA  +  xs)2  +  2  Mgr2  (A-13) 

The  elements  on  the  mass  matrix  can  be  re-written  by  introducing  non- 
dimensional  parameters.  These  are: 


2m/ 

m7 

(A- 14) 

Mx  =  MF  (1  +  y) 

(A- 15) 

Ms 

(A- 16) 

These  expressions  are  the  result  of  non-dimensionalizing  the  equations  of 
motion,  as  will  be  shown. 

Strain  Energy 

The  strain  energy  for  the  model  consists  of  two  parts,  that  due  to  the  wing 
bending  and  that  due  to  an  artificial  plunge  spring,  k.  The  plunge  spring  is 
present  so  that  the  effect  on  aircraft  stability  of  restraining  the  plunge  degree- 
of-freedom  can  be  studied.  The  strain  energy,  U,  is: 

u  =  2  [|  JT '  EI  (f '  f  dy]  +  |  kw”-  (A-17) 

The  wing  cantilever  natural  frequency,  u0,  is  introduced  into  the  strain  energy 
by  noting  that,  for  free  vibration  of  the  clamped  wing  only, 

J  (F)2/0  EI(f '  )2  dy  =  ^  E2  a>2  /o  mf2dy  (A-18) 


or 
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|  (E)2  (  El  (P '  )2  dy  =  ±  mu2  (E)2  (  f2  dy 


thus, 

U  =  2  {■-  mw02  JT  f2  dyjE2  +  —  kw2 


(A- 19) 


(A-20) 


Generalized  Forces 

The  aerodynamic  forces  acting  upon  the  aircraft  arise  from  the  distributed 
airloads  along  the  wings  and  the  load  from  the  canard  surfaces.  These  loads 
are  deformation  dependent.  Unsteady  aerodynamic  theory  was  used  to  develop 
expressions  for  the  distributed  wing  lift  and  moment.  These  expressions 
according  to  Reference  4,  are  first  written  in  terms  of  wing  deflection,  H,  pitch, 
a,  and  derivatives  a  —  dH/dy,  and  r  =  da/dy.  These  expressions  are: 

L(y,t)  =  nph 2  {H  +  VdcosA  +  Vtr  sinA  -  bad 

-  ba  Vr  sinA}  +  2?r  p\>  V  cosA  C(k)  {H  +  V  a  cosA 

+  Vtr  sinA  +  b  (— — a)  (a  +  VYsinA)}  (A-21) 

2 

My(y,t)  =  -jrpb3  { V  ( —  -  a)  a  cosA  +  —  V2tcosA  sinA 

2  2 

—  aH  —  Va  trsinA  +  b  (—  +  a2)  (d  +  Vr  sinA)} 

8 

+  2  7rpb2  VcosA  (~  +  a)  C(k)  {H  +  V  a  cosA 
& 

+  Vtr  sinA  +  b  (— ■  —  a)  (d  +  Vr  sinA)} 

mt 


(A-22) 


« 
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+  Vo  -  V  sinA  42.}  +  2  npb  (VWA)  C(k) 
dy 

{(x  +  y  sinA  +  b  (  —  -  a)  cosA)  0/(V  cosA) 

2 

-  h/(V  cosA)  -  w/(V  cosA)  +  0/cosA  -  tanA  (A-27) 

dy 

My(y,t)  =  -Jty?b3  {-Va0  4-  (V/2)  cos2A  0  +  aw  4-  ah 

-  (x  +  y  sinA)  a 6  4-  Va  -7^-  sinA 

dy 

4-  b  4-  a2)  0  cosA} 

O 

4-  2xpb2  V2  cos2A  C(k)  {0/cosA  -  tanA}  4-  a) 


I 
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l 


< 


i 


+  2rpb2V  cosA  +  a)  C(k)  {-h  -  w  +  (x  +  y  sinA) 

mt 

+  b  (-  -  a)  cosA)0}  (A- 28) 

mt 


Virtual  Work: 

If  the  aircraft  is  given  an  arbitrary  virtual  displacement,  consisting  of 
6h,  6w,  60,  the  virtual  work  done  by  the  airloads  on  the  wing  is  written  as 
follows: 

6  W  =  [/  {-L(x+y  sinA)  +  M  cosA}  dy]60  4-  If  fLdy]  6h 

+  [/  Ldy]  6w  (A-29) 

The  generalized  forces  to  be  used  in  the  equations  of  motion  are  found  from 
equation  A-29.  They  are: 


Q,  =  J  (-L  (x  +  y  sinA)  +  M  cosA)  dy 

(A-30) 

Qh  =  (  n-dy 

(A-3I) 

Q»  =  }'  Ldy 

(A-32) 

Substituting  the  expressions  for  L,  M,  and  f,  and  then  performing  the 
required  integrations,  the  expressions  for  the  generalized  forces  become: 

Qw  =  irpb2  {(x/  +  —■  /2  sinA  -  ba /  cosA)  0  -  — -  /  R 
J  o 

-  /w  +  V/0-  (VsinA)  R) 


I 
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Qw  =  7r/>b2  {(x/  +  —  /2  sinA  -  ba /  cosA)  6  —  —  /  E 
2  5 

-  /w  +  V/0-  (VsinA)  E) 


+  2xphV  cosA  C(k)  {(x/  +  -—/2  sinA  +  hi  (— —  a)  cosA)  9 

2  2 

-  ~  it  -  /w  +  l\9  -  (V  sinA)E} 
o 


Qh  =  npb~  l  {(-—x  +  —  /  sinA  -  —  ba  cosA)0 
5  45  5 


,  104  .  r-  2  ,  2  1  ..  .  ,  h  i 

+  T  ve  ~  7  v  smA  7) 


+  2npbl  (V  cosA)  C(k)  [(•— x  +  ■—  /  sinA 

5  45 

.i/I  ,2  k\  a  104  r-  2 
+  b(7-a)IoosA)9-  — E--w 

+  f  V9-(i  VsinA)E| 

Q,  =  —xpb'1  {(x2/  +  x/2  sinA  +  l3  sin2A)  -  (ba  cosA) 

o 


(x/  +  ~  /2  sinA)]  0  -  (“x/  +  -^7  /2  sinA)  E 
2  5  45 

-  (x/  +  ■£*  l~  sinA)  w  +  (x/  +  ~  /2  sinA)  V0 
&  £ 


-  (Vx  sinA  +  7-  /V  sin2A)  E} 
5 


-  2  Tr^bV  cosA  C(k)  { [(x2/  +  x/2  sinA  +  /3  sin2A) 

O 


(A-33) 


(A-34) 
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-  (Vx  sinA  +  —  V7  sin2A)  E} 


-  x/)b3  [(—  V7  cos3A  —  V/a  cosA)  0  +  (a/  cosA)  w 
2 


2  2 

+  ( —  a /  cosA)  E  —  (x/  cosA  4 - /2  sinA  cosA)  a0 

o  2 


+  (Va  sinA  cosA  (— ))  E  +  (b /  (—  +  a2)  cos2A)  0] 
5  8 


+  2ff/?b2  V  cosA  (—  +  a)  C(k)  {— (—  /cosA)  E  (A*35) 

2  5 


-  (/  cosA)  w  +  [(x/  +  —  /2  sinA)  cosA  +  b/  (—  -  a)  cos2A)  0} 

2  2 


The  equations  of  motion  for  the  airplane  are  now  written  as 


[Mij] 


w 

w 

Qw 

E 

+  lkijl  ' 

E 

—  , 

0 

0 

Q# 

(A-36) 


Where  k,j  =  0,  except  for  the  following  two  terms,  found  from  the  strain 
energy  expression,  U. 


.  _  d2U  _  , 

k"  '  "  k 


.  _  d2U  _  „  o  2  , 

koo  -  — —  -  2mw.‘J  f*  dy 


(A-37) 

(A-38) 


A  set  of  non-dimensional  coordinates,  f,  are  defined  by  the  following 
transformation. 
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►  < 

1  //  0  0 

w 

u 

~ 

0  l/l  o 

E 

a* 

0  0  1, 

9 

(A-39) 


Inverting  this  relationship  we  have: 


w 

1  0  0 

E 

— 

0  /  0 

*2 

9 

0  0  1, 

£3 

=  M 


(A-40) 


Equation  A-36  is  transformed  as  follows: 

[m]  [ip]  {£}  +  [ip]r  [k]  [ip]  {£}  =  [ip]T  {Q}  (A-41) 


Each  equation  (row)  is  then  divided  by  MT/2.  Equation  A-41  becomes: 


[My]  {{>  +  Iky)  {{}  =  (1/Mt/3)  MT  {Q}  (A-42) 


At  this  point  we  use  notation  (“)  above  a  variable  to  indicate  division  by  /,  e.g. 


7  =  z/l  JTS  =  xjl 


Non-dimensional  terms  appearing  in  the  equations  of  motion  are:  /i  =  wing 
mass/fuselage  mass,  f}  —  store  mass/wing  mass.  The  non-dimensional  elements 
of  the  mass  matrix  are: 


M„  =  1  +  1-^) 

(A-43) 

M-  =  l  <7+7>  + 

(A-44) 

M13  =  0 

(A-45) 
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s 

w 

n 

o 

(A-45) 

M22  =  —  (— )  +  (-^-) 

22  405  '  1  +/i  '  V  1-f/i ' 

(A-46) 

=  -(Ttr1 1 5 T  +  « sinA> 

(A-47) 

(7  +  sinA  +  7S) 

Ms3  ■  7T7  +  " (7  +  7  s5nA)2  +  ( 7+7*  (17  sin°‘A  + 

r‘  cos2A)  +  (^jj)  T2  +  (y^)  (A-48) 

[(7  +  -y  sinA)2  +  (7  +  “  sinA)  (*B  +  —  sinA)] 


+ 


2,3° 

£jt 

1+/* 


(7  +  sinA  +  JT6)2 


Mf(1  +/i)/2 


(A-49) 


The  non-zero  terms  of  the  non-dimensional  stiffness  matrix  are: 

kjj  k/Mx  —  utp  (A-50) 

where  wp  is  the  uncoupled  plunge  natural  frequency. 

k2o  =  w2  (2(m//Mx)  JT  f2  dy]  =  w2  (y^-)(~|)  (A-51) 

The  generalized  forces  can  be  divided  into  two  categories,  the  apparent  mass, 
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following  parameters,  we  can  write  these  apparent  mass  and  AIC  terms  in 
matrix  form. 


2  bl  =  S 


D  = 


pSCLa  Vn 

Mt 


a  — - 

9 


Mn 


Q  = 


2qnSC1 

Mt 


7  =  (*  +  “  sinA) 


V  cosA  =  V„ 


The  elements  of  the  apparent  mass  matrix  are: 


11 

s 

(A-52) 

-  2  - 
M12  ~  J  » 

(A-53) 

Mi3  =  "7*  (7  +  “  E  cosA) 

(A-54) 

w  -  104  - 

M"  "  405  '' 

(A-55) 

M =  -p[y  +  ~  sinA  +  ~  C  cosA) 

(A-56) 

M33  =  ~li  (r  +  sin2A  +  yF  cosA  +  cos2A) 

(A-57) 

M21  =  M12 

(A-58) 

^31  =  ^13 

(A-59) 
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M30  —  M03 

The  elements  of  the  aerodynamic  influence  coefficient  matrix  are 
K;j  +  io;b;j.  The  elements  of  this  matrix  are  listed  below: 

Kn  -  K21  =  K31  =  0 

Kjo  =  QC(k)/tanA 

Kj3  =  “QC(k)/cosA 

Koo  —  “  QC(k)  tanA 

mt 

K03  =  ~~  QC(k)/cosA 
5 

K32  =  -QC(k)  tanA  (5T  +  sinA) 

o 

K33  =  QC(k)y/cosA 
b„  =  DC(k) 

b12  =  D  ("rE  tanA  +  -7  C(k)) 

«  O 

bi3  =  [7  &  +  C(k)  (7  +  E  cosA)] 

b21  =  t  DC(k) 

5 

boo  =  D[—  E  tanA  +  (  — )  C(k)l 

b23  =  “D[t  — ~  +  c(k)  T  (7  +  -r  SinA  +  E  cosA)] 

5  cosA  5  9 


(A-60) 
given  as 

(A-61) 
(A- 62) 
(A-63) 

(A- 64) 

(A-65) 

(A-66) 

(A-67) 

(A-68) 

(A-69) 
(A-70) 
(A- 7 1 ) 
(A-72) 

(A-73) 
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bji  =  -DC(k)  7 

b32  =  -D[-^-  E  tanA  (7  4-  ~  sinA)  + 

“  C(k)  (7  +  ■“  sinA)  +  j  E2  sinA] 

b33  =  D[y  (7/cosA)  +  C(k)  (T2  +  “  siu2A  +  E  7  cosA) 

— —  E2  sin2A] 

2  ' 


(A-74) 

(A-75) 


(A-76) 
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APPENDIX  B 
Solution  Methods 


Two  methods  were  used  to  solve  the  flutter  problem.  One  method,  called 
the  P-K  method,  assumes  a  system  response  both  at  flutter  and  before  and 
after  flutter  of  the  form  ept,  where  p  is  a  complex  number.  For  simple 
harmonic  motion  at  flutter  the  aerodynamic  influence  coefficients,  which  are 
derived  from  the  generalized  aerodynamic  forces,  are  valid.  With  the  P-K 
method,  the  eigenvalue  equation  is  written  as: 

|p2|M|  +  (K|  -  <UA1|  Ui)  =  {0}  (B-l) 

The  aerodynamic  influence  coefficients  can  be  divided  into  real  and  imaginary 
parts.  They  are  written  as: 

qnIAl  =  qn  IQrJ  +  p(~)  [Qil  (b-2) 

tv 

where  [  QR  j  and  [  Qj  ]  are  the  real  and  imaginary  parts  of  [A],  respectively. 
Substituting  Eqn.  (B-2)  into  Eqn.  (B-l),  we  have: 

|p2|M|  +  |K|  -  q„IQR|  -  p  (— )  [Q,||  «,}  =  {0}  (B-3) 

C v 

or 


Ip2[m]  +  P(B]+IC|]  {$}  =  (0} 


(B-4) 


< 
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where 

|B]  =  -(-^-1  IQil  =  JfiV,  (£)  [Qi] 

and 


(B-5) 


(B-6) 


(B-7) 


The  method  begins  with  the  choice  of  an  initial  value  of  reduced  frequency 
and  a  fixed  value  of  velocity.  This  initial  reduced  frequency  is  a  starting  point 
for  calculating  the  aerodynamic  influence  coefficients.  A  value  of  p  is  found 
from  Eqn.  (B-5).  Recall  that  p  is  a  complex  number. 


P  “  Pr  +  »Pi 


(B-8) 


From  the  imaginary  part  of  p,  a  new  value  for  reduced  frequency  can  be 
calculated.  This  new  value  of  reduced  frequency  generally  will  not  match  the 
initial  guess  for  reduced  frequency.  Therefore,  a  “lining  up”  process  of  the 
input  values  k  with  the  computed  values  of  k  is  necessary.  This  lining  up 
process  continues  until  the  input  value  of  k  matches  the  output  value  of  k  for  a 
particular  frequency.  This  process  is  then  repeated  for  each  modal  frequency 
and  for  each  velocity  increment. 

The  other  method  used  to  solve  the  flutter  problem  is  similar  to  the  P-K 
method.  Instead  of  using  the  unsteady  Theodorsen  aerodynamic  formulation, 
as  was  used  in  the  previous  method,  quasi-steady  aerodynamics  are  used.  In 


quasi-steady  aerodynamics,  motion  occurs  very  slowly  such  that  reduced 
frequency  can  be  taken  as  zero. 

With  this  assumption,  the  eigenvalue  equation  is  greatly  simplified.  Since 
reduced  frequency  is  zero,  the  [  QR  ]  and  [  Qj  ]  matricies  are  independent  of 
reduced  frequency.  The  resulting  eigenvalue  equation  is  written  as: 

[p2[M|  +  p[B|  +  [C]|  {{;}  =  {0}  (B-0) 

where 

IB]  =  -j  pbV„  [Q,l  (B-10) 

|C]  =  IK)  -  q„[QB]  (B-ll) 

The  f  Qr]  and  [  Qj]  matricies  are  the  real  and  imaginary  parts  of  [A]  with  the 
reduced  frequency  dependence  removed. 

Since  Eqn.  (B-10)  is  independent  of  reduced  frequency,  no  frequency 
iteration  is  necessary,  values  for  p  can  be  calculated  directly  for  a  specific 
velocity. 


